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Abstract 

Pulmonary surfactant is required for lung function at birth and throughout postnatal life. Defects in the surfactant system 
are associated with common pulmonary disorders including neonatal respiratory distress syndrome and acute respiratory 
distress syndrome in children and adults. Lipogenesis is essential for the synthesis of pulmonary surfactant by type II 
epithelial cells lining the alveoli. This study sought to identify the role of pulmonary epithelial SREBP, a transcriptional 
regulator of cellular lipid homeostasis, during a critical time period of perinatal lung maturation in the mouse. Genome wide 
mRNA expression profiling of lung tissue from transgenic mice with epithelial-specific deletions of Scap (Scap^^'^, resulting 
in inactivation of SREBP signaling) or Insigl and Insig2 {lnsigl/2^'^, resulting in activation of SREBP signaling) was assessed. 
Differentially expressed genes responding to SREBP perturbations were identified and subjected to functional enrichment 
analysis, pathway mapping and literature mining to predict upstream regulators and transcriptional networks regulating 
surfactant lipid homeostasis. Through comprehensive data analysis and integration, time dependent effects of epithelial 
SCAP/INSIG/SREBP deletion and defined SCAP/INSIG/SREBP-associated genes, bioprocesses and downstream pathways 
were identified. SREBP signaling influences epithelial development, cell death and cell proliferation at El 7.5, while primarily 
influencing surfactant physiology, lipid/sterol synthesis, and phospholipid transport after birth. SREBP signaling integrated 
with the Wnt/p-catenin and glucocorticoid receptor signaling pathways during perinatal lung maturation. SREBP regulates 
perinatal lung lipogenesis and maturation through multiple mechanisms by interactions with distinct sets of regulatory 
partners. 
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Introduction 

Pulmonary surfactant forms mono- and multilayer films that 
line the alveolar surface, reducing the surface tension created at 
die air-liquid interface to ensure proper inflation and function of 
the lung. Qualitative and quantitative defects in the surfactant 
system are associated with common pulmonary disorders includ- 
ing neonatal respiratory distress syndrome (RDS), bronchopulmo- 
nary dysplasia (BPD) and acute respiratory distress syndrome 
(ARDS) in children and adults. Mutations in genes critical for 
surfactant production or function, including surfactant proteins A 
{SFTPA), B {SFTPB) and C {SFTPQ, a.nd ABCA3 are associated 
with chronic lung disease in humans and mice [1,2]. 

Surfactant is a protein -lip id complex comprised of 90% lipid 
and 10% surfactant-associated proteins. Phospholipids represent 
70-80% of the surfactant by weight, of which 45% is represented 
by disaturated phosphatidylcholine (SatPC) [3]. Transcriptional 
induction of pathways regialating lipid synthesis and surfactant 
proteins occurs during the later stages of perinatal lung maturation 
and is required for respiratory function at birth; thus, alveolar lipid 



and surfactant pool size requires precise regulation to ensure 
proper adaptation to the transition to air breathing. 

Sterol regulatory element -binding proteins (SREBPs) are master 
transcriptional regulators of cellular lipid homeostasis active in 
multiple tissues. Three distinct SREBP isoforms, SREBP- la, -Ic 
and -2, are encoded by two genes [Srebfl and Srebf2). In general, 
SREBP- la and SREBP- Ic regulate fatty acid synthesis while 
SREBP-2 regulates cholesterol synthesis [4,5]. When intracellular 
cholesterol levels are abundant, SREBP proteins are held in an 
inactive form in the endoplasmic reticulum (ER) complexed to the 
sterol-sensing protein SCAP (SREBP cleavage-activating protein) 
and the ER-resident proteins INSIGl and INSIG2. In sterol- 
depleted cells, INSIGl/2 dissociate from SCAP and the SREBP- 
SCAP complex is transported to the Golgi where the active form 
of SREBP is generated and released from the membrane by two 
sequential proteolytic cleavages by Site-1 and Site-2 proteases 
[6,7]. The active form of SREBP enters the nucleus where it 
regulates transcription of target genes. 

Lipogenesis is critical for the synthesis of pulmonary surfactant 
by type II cells. Phospholipid synthesis in type II cells was 
markedly increased shortiy before birth in preparation of air 
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breath. SREBP- Ic is expressed in the developing lung and its 
expression increases during late gestation, in association with 
increased surfactant lipid synthesis [8]. SREBP- Ic activates 
PCYTIA [9-11], the rate-limiting enzyme for phosphatidylcho- 
line (PC) synthesis, thereby increasing PC synthesis in type II cells. 
Our recent studies showed that deletion of Scap in the adult mouse 
lung inhibited SREBP activity in respiratory epithelial cells, 
resulting in altered pulmonary lipid homeostasis associated with 
compensatory lipid synthesis and accumulation in lipofibroblasts 
[12]. Conversely, deletion of the Insigl and Imig'2 genes (Insigl 12^^ 
^) acti\'atc(l SREBP-mediated lipogenesis in rc'sj)irat()ry c'pithelial 
cells resulting in lipotoxicity-related lung inflammation and tissue 
remodeling in adult mice [13]. Transcriptional network modeling 
revealed that SREBP is an important regulatory hub in a lung lipid 
transcriptional network [14]. Collectively, data from both 
bioinformatic and mouse models demonstrated important roles 
for the SREBP signaling in lung lipid homeostasis. 

This study was designed to further define the dynamic roles of 
SREBP in perinatal lung development using integrative approach- 
es. Through comprehensive data analysis and integration, we 
identified early vs. late effects of epithehal SCAP/INSIG/SREBP 
signaling and defined SCAP/INSIG/SREBP-associated genes, 
bioprocesses and downstream pathways. We show that SREBP 
signaling influences epithelial development, cell death and cell 
proliferation at El 7.5, while primarily influencing surfactant 
physiology, lipid/sterol synthesis, and phospholipid transport at 
PNl. In addition to regulation of the lipogenesis program that was 
dominant during the later stages of lung development, our study 
demonstrates clear cross-talk between the SCAP/INSIG/SREBP 
pathway and two other signaling pathways important for lung 
development and maturation, specifically Wnt/ |3-catenin and 
glucocorticoid receptor signaling. 

Materials and Methods 

This study was carried out in accordance with the recommen- 
dations in the Guide for the Care and Use of Laboratory Animals 
of the National Institutes of Health. All study protocols were 
approved by the Institutional Animal Care and Use Committee of 
the Cincinnati Children's Hospital Rcs(^ar(ii Foundation (Permit 
Number: 2C12114). All microarray data are MIAME compliant 
and submitted to GEO (http:/ /www.ncbi.nbn.nih.gov/geo/). 

Transgenic Mice 

ScafP""'^"' mice containing a loxP-flanked neomycin resistance 
cassette located 3 kb upstream of exon 1 and an additional loxP 
site in intron 1 were obtained from Jackson Laboratory' (Bar 
Harbor, ME). Epithehal-specific, conditional Scap knockout mice, 
termed Scap'^'^ , were generated crossing Sca^"'^' mice with: 1) 
mice that express rtTA protein in the respiratory epithelium (SP-C- 
rtTA~'% and 2) tetO-Cre recombinase mice {(tetO)yCMV-Cre'^'% 
This colony was maintained in a mixed FVB/N and B6/129S6 
background. 

Sca^"'^'" Uttermates lacking either rtTA or Cre genes served as 
controls. Dams were administered doxycycline chow (625 mg/kg, 
Harian Teklad, Madison, WI) from E6.5 to El 2.5 to delete Scap 
from respiratory epithelial cells. Importantly, doxycycline was 
withdrawn at El 2.5 to circumvent possible rtTA/Cre-mediated 
toxicity effects. Lung tissue was collected from pups at El 7.5, 
E18.5 and PNl (day of birth). 

Insigt^"'^" mice have LoxP sites flanking exon 1 and Insig2 
knockout mice (Insig2 ^ ) were generated by replacing exons 2 
and 3 with a neo cassette. Insigl^''-^' and Insig2 ' mice were 
obatined from The Jackson Laboratory (Bar Harbor, ME). 



Doxycycline-exposed (E6.5-12.5) SP-C-rtTJt / {ktO)jCMV'Crt ' 
'^llnsigJ^''^' IInsig2~'~ mice are termed Insigl /2^'^ mice {Insigl 
was conditionally deleted in the Insig2 ' background). Insigl^'"'' 
fi"IInsig2~'~ Uttermates lacking either SP-C-rtTA or (tetO)yCMV- 
Cre served as controls. Genotypes were identified by PCR from 
genomic tail DNA as described previously [15,16]. 

Mice were maintained in a pathogen free environment under 
the supervision of the Institutional Animal Care and Use 
Committee of Cincinnati Children's Hospital Research Founda- 
tion. Mice were kept on a 14-10 light-dark cycle and were allowed 
access to food and water ad libitum. 

Tissue Collection 

E17.5 and E18.5 mice were delivered by Cesarean section with 
anesthesia followed by exsanguination. PNl mice were taken on 
the day of birth and euthanized with pentobarbital followed by 
exsanguination. Lungs were flash frozen for analysis. 

Cell Transfection 

MLE15 cells were cultured in HITES media [17]. For study, 
5 X 1 0'^ cells were plated in 6-well plates the day before 
transfection. Cells were transfected with 30 pmol/weU of iiScap 
(Silencer Select, Life Technologies Corporation) using Lipofecta- 
mine RNAiMAX Reagents (Life Technologies Corporation) 
according to the manufacturers protocol. Cells were harvested 
48 hr post-transfection for RNA isolation and qPCR analysis. 

RNA Isolation and qPCR 

Total RNA was isolated from whole lung homogenates using 
the RNeasy Mini Kit (Qiagen, Valencia, CA) and reverse 
transcribed using the iScript kit (BioRad, Hercules, CA). miRNA 
was isolated using the miRNeasy Mini kit with on-column DNase 
digestion (Qiagen, Valen[:ia, CA) and reverse transcribed using the 
Taqman MicroRNA Reverse Transcription Kit (Life Technolo- 
gies, Grand Island, NY). Quantitative PCR was performed on a 
StepOnePlus Real Time PCR System (Applied Biosystems, Foster 
City, CA). Data for RNA analysis were normalized to 18 s and 
data for miRNA analysis were normalized to snoRNA202. 

Saturated Phosphatidylcholine (Sat PC) and Total 
Phospholipid IVleasurement 

Total lipids were extracted from fetal lung tissue using 
chloroform and methanol according to the Bhgh and Dyer 
method [18]. Sat PC was then isolated using osmium tetroxide and 
neutral alumina by the Mason, et. al. method [19]. Isolated SatPC 

and total phospholipid were then quantitated by phosphorus 
measurement as previously described [20]. 

Cholesterol and Triglyceride Quantification 

Total cholesterol, cholesterol esters and triglycerides were 
determined using the Cholesterol/Cholesterol Ester Quantifica- 
tion kit (EMD Chemicals, Inc) and the Triglyceride Quantification 
kit (Bio Vision) according to the manufacturer's protocol. 

Microarray Analysis 

Total lung RNA was isolated from Scap'^^^, Insigl 1 2^'^, and 
respective control Uttermates at E17.5, E18.5 and PNl using 
RNeasy Mini Kit (Qiagen, Valencia, CA). The cRNA was then 
hybridized to the Mouse Gene 1.0 ST Array (Aflymetrix, Santa 
Clara, CA) according to the manufacturer's protocol. Three 
biological replicates were used for each condition. The RNA 
quality and quantity assessment, probe preparation, labeling, 
hybridization, and image scan were carried out in the Cincinnati 
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Children's Hospital Medical Center Affymetrix Core using 
standard procedure. Data were further analyzed using Gene- 
Spring GX 1 1 .5 (Agilent Technologies, Santa Clara, CA). 
DifTerentially expressed genes in Scap'^'^ vs. control mice at 
different developmental time points were identified using two-way 
ANOVA (genotype and time dependc-nt gene changes). Genes 
altered in vs. control were identified using Student's t- 

test with the threshold of corrected PsO.Ol and fold change Si. 5. 
Benjamini-Hochberg false discovery rate (FDR) was applied for 
multiple testing correction. The probe sets were pre-filtered by 
their expression range within 35-100 percentiles for at least 66% 
of samples to remove probes with extremely low expression signals 
in both mutant and control samples. The complete dataset has 
been submitted to the Gene Expression Omnibus (GEO) database 
with the accession number of GSE53397. 

miRNA Analysis 

miRNA expression profile was determined by using GeneChip 
miRNA 3.0 Array (www.aHymetrix.com). This array platform 
provides complete coverage of miRBase vl7, snoRNA and 

scaRNA probes on the same array, and covers 153 organisms 
including human, mouse and rat. Lung samples from Scap^^^ and 
control littermates at £18. 5 (n = 3 per genotype) were used to 
isolate RNA using miRNeasy kit with on column DNase digestion 
(Qiagen, Valencia, CA). The Expression Console software 
(Afi"ym(;trix, Inc.) was used for data summarization, normalization, 
and quality control. The miRNAs with P<0.01 and fold changes 
> 1 . 3 was defined as differentially c-xpressed u.sing BRB-Array- 
Tools Version: 4.2.1 (http://Unus.nci.nih.gov/BRB-ArrayTools. 
html). miRNA target prediction and filtering were performed 
using miRNA Target Filter tool in IPA. miRNA and mRNA 
targets showing inverse expression patterns in response to Scap 
deletion at E18.5 were selected for further analysis (i.e. miRNA up, 
mRNA down OR miRNA down, mRNA up). 

DiflFerentially expressed genes at difi^erent lung maturation 
stages (i.e. El 7.5, E18.5 and PNl) were subjected to Self- 
organizing Map (SOM) clustering, gene ontology analysis, 
promoter analysis, functional classification and pathway analysis 
as described previously [21,22]. Gene Ontology and Pathway 
analysis were performed using publicly available web-based tool 
DAVID (database for annotation, visualization, and integrated 
discovery), ToppGene Suite (http://toppgene.cchmc.org/) and 
Ingenuity Pathway Analysis (IPA, Ingenuity). 

Upstream regulator analysis uses Ingenuity software based on 
prior knowledge of expected effects between transcriptional 
regulators and their target genes pre-compUed in the Ingenuity 
Knowledge Base. The analysis examines how many known targets 
of each transcriptional regulator are present in the user's dataset, 
and also compares their direction of change in the experimental 
condition to what is expected from the literature. If the observed 
direction of change is positively correlated with the literature, then 
we can predict the transcriptional regulator is at activated state. 
For each potential transcriptional regulator, two statistical 
measures (overlap p-value and activation z-score) are computed. 
The p-value is calculated to rank upstream regulators based on 
significant overlap of the differentially expressed genes to the 
known targets regulated by a transcriptional regulator. The z-score 
is used to infer likely activation states of upstream regulators. 

Promoter sequences (2 kb upstream of the transcriptional start 
site) were retrieved using Biomart (www.ensembl.org/biomart/ 
martview). Multiple promoter analytic tools (Genomatix Matin- 
spector. Clover, CisView, and Paint4) were used to identify 
enriched transcription factor binding sites in the promoter regions 
of genes in the cluster of interest. Cytoscape (http:/ /cytoscapeweb. 



cytoscape.org/) was used to map associations between transcrip- 
tion factors matrix families and predicted target genes. PBGE 
database (http://research.cchmc.org/pbge/) was used to match 
the transcription factors matrix families to the transcription factors 
expressed in lung. 

Results 

The Effects of Epithelial Scap Deletion on Perinatal Lung 
Lipid Content 

Mice with conditional deletion of the Scap {Scap'^^^) or Insigl 
{Insigl /2^^^) genes in pulmonary epithelium were generated as 
previously reported [12,13]. Scap and Insigl gene deletion 
efficiency was assessed by qPCR and exon deletion was confirmed 
by exon microarray, demonstrating a marked decrease in Scap 
(Figure lA) and Insigl (Figur(- IB) mRNA in our mouse models. 
We have previously demonstrated that epithelial SCAP/ SREBP 
signaling regulates lipid homeostasis in the adult lung [12,13]. To 
evaluate the effects of epitheUal Scap gene deletion on perinatal 
lung lipid homeostasis, we measured lung SatPC, total cholesterol 
and triglyceride levels in lungs at E17.5, E18.5 and PNl. SatPC 
was significandy decreased at E18.5 and PNl in Scap^'^ mice vs. 
control littermates, but not at E17.5 (Fig. 2A). Changes in SatPC 
content correlated with the deletion efficiency oiScap in this model 
(Figure lA). Both total cholesterol and triglyceride content were 
not significantly altered by group mean comparison bc-twecn 
Scap and control mice at any of the time points examined 
(Fig. 2C and E). Considering the sample variation related to the 
incomplete deletion efficiency of the SP-C/btO-Cre transgenic 
system, we examined the correlation between Srehf-lc expression 
levels and lung lipid content at PNl. Results from this analysis 
showed that Srebf-lc levels positively correlated with SatPC and 
total cholesterol content, with correlation coefficients of 0.63 and 
0.67, respectively; however, there was no strong correlation 
between Srebf-lc and triglyceride levels (—0.27) (Fig. 2B, D and 
F). These results indicated that, due to the variation of deletion 
efficiency and the possible compensatory effect of other cell types 
in the lung to the epithelial deletion of Scap in vivo, the t-test 
comparison of group mean may not be sensitive enough to reveal 
the direct effect of Scap deletion in epithelial cells on cholesterol 
levels. These data demonstrate that epithehal SCAP/SREBP 
signaling positively regulates SatPC content during perinatal lung 
maturation. 

The Functional Genomic Responses to SREBP Deletion/ 
Activation during Perinatal Lung Maturation 

To systematically determine the genomic responses of SREBP 
deletion/activation in the perinatal lung, genome wide mRNA 
expression profiles of Scap'^^^ mouse lungs were performed at 
E17.5, E18.5 and PNl; Insigl / 2^^" mouse lungs were analyzed at 
El 7.5. Genes differentially expressed in Scap"^'^ mice, Insigl 12"^ 
mice and their corresponding control littermates in at least one of 
the conditions were identified by ANOVA. This analysis identified 
184 genes that were up-regulated and 1079 genes down-regulated 
in Scap^'^ mice, suggesting that SREBP works primarily as a 
transcriptional activator in the perinatal lung (Figure 3A and 
Table SI). Among the genes that increased in any of the three time 
points in Scap"^'"^ lungs, "protein targeting to ER/membrane", 
"RNA catabolic/metaboUc process" were the most enriched 
functional categories (Figure 3B). Among the genes that were 
suppressed in Scap'^^^ lungs, "regulation of transcription", 
"epithelium development", "vesicle mediated transport/sorting", 
"lung development", "lipid biosynthetic/metabohc process", and 
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Figure 1. Relative expression of Scap and Insigl in Scap^^^ and mice. Scap (A) and Insigl (B) mRNA were measured by mRNA 

microarray analysis and qPCR. Total lung RNA was isolated from Scap''^^, and littermate controls at E17.5, E18.5 and PN1. Data were 
normalized to 18S and presented as mean ± S.E.; n = 2-3/group. 
doi:1 0.1 371 /journal.pone.0091 376.g001 



"regulation of cell death/cell cycle" were most enriched functional 
classes (Figure 3C). 

We further assessed the functional enrichment of genes altered 
at El 7.5 (early effect) and PNl (later effect) separately to identify 
the dynamic effects of SREBP deletion/activation on perinatal 
lung development and maturation. Genes involved in "Cellular 
movement", "Cardiovascular System Development and Func- 
tion", "Respiratory System Development and Function", "Cellu- 
lar Assembly and Organization", "Cell Morphology", "Cell 
Cycle" were the most highly enriched bioprocesses at E17.5, 
decreasing with time in response to disruption of SREBP signaling. 
In contrast, genes involved in "Lipid/Sterol Synthesis/Metabo- 
lism", "Angiogenesis", "Epithelial development/ differentiation", 
were mostly highly enriched at PNl (Table 1). Canonical 
pathways significantiy affected at El 7.5 by Scap deletion include 
"Oxidative Phosphorylation", "Mitochondrial Dysfunction", 
"Wnt/ P-catenin Signaling", and "Glucocorticoid Receptor Sig- 
naling"; at the PNl time point, pathways that were enriched 
included "Steroid Biosynthesis", "Ephrin Receptor Signaling", 
"Integrin Signaling" and "LXR/RXR Activation" (Table 2). 
Many of the known targets regulated by SREBP signaling were 
decreased in Scap'^'^ mice and/or induced in Insigl 1 2^'^ mice 
(Table S2), including Scdl, Scd2, Fasn, Ldlr, Hmgcsl, Sqle and Lss, 
demonstrating an inverse relationship between SREBP suppres- 
sion and activation in Scap'^'^ and Insigl 12^^^^ mice, respectively. 
In addition to the known effects of SREBP signaling on lipid 
metabolism that was observed predominandy at PNl, genes 
{Sftpal, Sftpb, AbcaS, Ctnnlb and NrScl) and pathways (Wnt/p- 
catenin Signaling) important in fetal lung development, matura- 
tion and function were decreased in Scap'^^^ mice at El 7.5. qPCR 
validation of a subset of these microarray data are shown in 



Figure 3D and 3E. Collectively, these data demonstrate dynamic 
roles of epithelial SREBP signaling during perinatal lung 
development. 

In order to understand the regulatory mechanisms and identify 
potential cofactors of the SREBP signaling pathway in regulating 
lung lipid homeostasis, genes differentially expressed in Scap'^'^ 
mice vs. control from two-way ANOVA analysis was clustered into 
8 different sub-groups based on expression pattern similarity using 
the SOM clustering method (Table SI, Figure 4A). Genes within 
each cluster were subjected to functional classification using 
Toppgene suite (http:/ /toppgene. cchmc.org/). Gene set enrich- 
ment analysis identified two "lipid" enriched clusters, 5 and 7. 
Cluster 5 genes were functionally enriched in molecular transport 
(Slc25al, Slcllal, AcslSj, protein trafficking {Arf6, Scamp2, Rab27d), 
lipid/phospholipids biosynthesis and metabolism (Pcytla, Lpcatl, 
Ljpla2, Scdl, Aacs, Stxl2) (Figure 4B). Genes within this cluster were 
moderately induced before birth in control mice but decreased 
within the same time frame in Scap^^^ mice (Figure 4A). Cluster 7 
genes were induced with age in control mice but maintained at a 
low level from E17.5 to PNl in Scap'^^^ mice. Genes within this 
cluster were functionally involved in cholesterol/steroid synthesis 
and in response to lipid [Hmgcsl, Mil, Sqle, Fdps, Msmol, Cdipt, 
Insigl, DgatT) (Figure 4C). An important functional diflFerence 
between the genes regulated at El 7.5 compared to PNl (in 
addition to regulation of lipid homeostasis) was the finding that a 
number of genes in cluster 7 play important roles during lung 
morphogenesis, maturation and surfactant function including 
Foxa2, FgflS, Sftpal, Sftpb and Abca3. 

To understand the role of SCAP/INSIG/SREBP signaling m 
the regulation of lung lipid homeostasis at a network level we 
sought to identify potential upstream regulators and downstream 
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Figure 2. Altered lung surfactant SatPC content in Scap^'^ mice. SatPC (A), total cholesterol (C) and triglycerides (E) were measured using 
whole lung samples from control (open bars) and Scap"^^^ mice (closed bars) at El 7.5, El 8.5 and PN1. *p<0.05. n = 5-8/group. Correlation coefficients 
between Srebf-lc mRNA and SatPC (B), cholesterol (D) and triglyceride levels (F) were measured at age PN1, n = 5-8/group. 
doi:1 0.1 371/journal.pone.0091 376.g002 



targets using promoter and literature mining approaches. Proxi- 
mal promoter sequences of genes in cluster 5 and 7 (2 kb upstream 
of each gene) were extracted and promoter analysis was performed 
to identify significantly over-represented cis-elements and corre- 
sponding transcription factors that are expressed in lung. A 
transcriptional network was constructed using the top 10 ranked 
transcription factor binding sites and their potential target genes 
(genes that contain positive binding sites in promoter region) in 



Cytoscape (http:/ /www.cytoscape.org/). As shown in Figure 5, 
transcription factors in the Lef, Stat, Cebp, Ebox/Creh, Fkhd, SPl, Pax 
and JVrfl families are important regulatory hubs in the network 
and likely to be partners of SREBP that regulate perinatal lung 
lipid homeostasis. Literature mining was performed using IPA 
knowledge base (https://analysis.ingenuity.com/). Potential up- 
stream regulators were predicted using genes from cluster 5 and 7 
respectively. The biological relationship of predicted core regula- 
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Figure 3. Genes and biological processes affected by Scap deletion in the developing lung. A) Differentially expressed genes in Scap 
lungs were identified by mRNA microarray analysis. Heatmap analysis shows the two-dimensional clustering of 1079 mRNAs that were significantly 
altered in Scap'^'^ mice compared with controls. The intensity in the red and green color ranges indicates increased versus decreased mRNAs, 
respectively. Each row represents a single mRNA, and each column represents an average fold change between knockout and control mice at the 
specified gestational age. B and C) Biological processes significantly influenced by epithelial Scap deletion were identified using ToppGene Suite 
(http://toppgene.cchmc.org). Enriched bioprocesses in genes induced (B) or suppressed (C) in Scap''^^ lungs are represented as red bars and green 
bars respectively. Statistical significance of each bioprocess was presented using negative log2 transformation of P value. D and E) qPCR analysis 
confirmed RNA microarray data demonstrating altered expression of genes involved in lipid biosynthesis/transport (D) and in lung development and 
maturation (E) in Scap''^^ mice (filled bars) compared with control (open bars), n = 3/group, *p<0.05. 
doi:1 0.1 371 /journal.pone.0091 376.g003 



Table 1. Enriched functions altered by Scap deletion In perinatal lung. 





El 7.5 Scaff'^'' 




El 8.5 & PN1 Scap^''' 




Functions 


p value 


Functions 


p value 


Cellular Movement 


2.79E-27 


Synthesis of sterol 


1.29E-14 


Cardiovascular System Development and Function 


8.96E-24 


Synthesis of lipid 


3.16E-12 


Cellular Assembly and Organization 


2.27E-23 


Synthesis/metaboiism of cholesterol 


i.oiE-n 


Cellular Function and Maintenance 


2.27E-23 


Regulation of transcription 


3.90E-1 1 


Tissue Morphology 


771E-23 


Proliferation of cells 


1 .28E-08 


Cellular Growth and Proliferation 


1.33E-22 


Cell death 


1.83E-08 


Cell Death 


1.49E-15 


Differentiation 


9.57E-07 


Cell Morphology 


6.04E-15 


Angiogenesis 


1.13E-06 


Respiratory System Development and Function 


1.83E-13 


Development of epithelial tissue 


3.57E-06 


Cell Cycle 


1.52E-10 


Fatty acid metabolism 


4.09E-06 
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Table 2. Enriched canonical pathways altered by Scap deletion in perinatal lung 





El 7.5 Scaf/^^ 




PN1 Scaff'^'' 




Pathways 


p value 


Pathways 


p value 


Oxidative Phosphorylation 


3.16E-14 


Biosynthesis of Steroids 


5.89E-08 


Mitochondrial Dysfunction 


2.00E-13 


Ephrin Receptor Signaling 


3.31 E-06 


Axonal Guidance Signaling 


6.31 E-11 


Integrin Signaling 


1 .OOE-05 


Ubiquinone Biosynthesis 


3.98E-09 


CXCR4 Signaling 


8.51 E-05 


Wnt/p-catenin Signaling 


1.91E-08 


Wnt/p-catenin Signaling 


9.33E-05 


Paxillin Signaling 


2.29E-08 


IL-8 Signaling 


1 .38E-04 


PTEN Signaling 


1 .78E-07 


ILK Signaling 


1 .07E-04 


Phospholipase C Signaling 


2.40E-06 


Insulin Receptor Signaling 


3.80E-04 


Notch Signaling 


1.23E-05 


mTOR Signaling 


4.68E-04 


Glucocorticoid Receptor Signaling 


7.59E-04 


LXR/RXR Activation 


8.32E-04 
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tor.s and their associated targets (based on literature) are illustrated 
in Figure 4B&C. As expected, SCAP/INSIG/SREBP were 
predicted as core regulators for both cluster 5 and 7 genes, but 
other core regulators in the network were distinct, indicating 
multiple gene-specific regulatory mechanisms are involved. 

Identification of Potential Upstream Regulators of SREBP 
Signaling in Perinatal Lung Development 

Differentially expressed genes in Scap^'^ mice were subjected to 
an upstream regulator analysis to identify potential mediators that 
drive the overall genomic responses caused by deletion or 
activation of SREBP signaling. From this analysis, we identified 
potential regulators of SREBP signaling at different time points 
during lung maturation. Table 3 lists the top 10 candidate 
regulators that are most likely to be involved in the SREBP 
signaling pathway at E17.5, E18.5 and PNl, ranked by p-value. A 
group of transcriptional factors (TFs) including Foxo3, IrfJ, E2f4, 
Ctnnbl, Pou4fl, GUI, and Gata3 (Gata6 may share similar group of 
downstream targets) were predicted to regulate transcriptional 
responses to both the deletion oi Insig and Scap at El 7.5, but not 
at El 8.5 or PNl. These transcription factors are known to 
positively regulate cell morphogenesis and proliferation, and 
negatively regulate cellular metabolic and biosynthetic processes, 
and programmed ceU death. At E18.5 and PNl, SREBFl/2, 
Foxo4, Nrlh2/3, Ppara and Pparg were predicted to be important 
regulators of the SREBP signaling. The activation status, either 
activated or inhibited, of the predicted upstream regulators in 
Scap'^^'^ and Insigl/2^^^ mice are listed in Table 4. The finding 
that Srebfl and 2 are two of the most inhibited TFs (associated with 
the lowest z-scores) in Scap^'^ mice, and the most activated TFs in 
Insigll2'^"^ mice, demonstrates the efficacy of this method to 
predict novel partners, co-factors and signaling molecules influ- 
enced by SREBP signaling in the perinatal lung. In addition to 
Srehfl and 2, the activation status of Pparg, Mrlh, Myc and Cdkn2A 
were inversely predicted in Scap'^^'^ and lnsigl/2^'^ mice, 
indicating that these factors are likely to be proximate downstream 
regulators of the SCAP/INSIG/SREBP cascade. The activation 
status of Cebpa, Ahr, Irfl, Mrfl and Hnf4a were predicted to be in the 
same direction for Scap'^'^ and lnsigll2'^"^ mice, suggesting that 
these factors are likely influencing or co-regulating SCAP/ SREBP 
signaling by more than one mechanism. 



Suppression of Wnt/|3-catenin Signaling by Scap Deletion 
at El 7.5 

From the microarray analysis, a number of Wnt/p-catenin 
signaling components, including ligands (Wnt4/5a/7d), receptors 
{Fzdl/2/4, Lrp 1/5/6, Cdhl/2/4/5) and effectors [CAnnbl) were 
decreased in Scap'^'^ mice, particularly at El 7.5 (Figure 6). In 
addition, several factors known to regulate Wnt-induced tran- 
scription, including NEMO-like kinase (jVft), Sox5/ 11 and retinoic 
acid receptor [Rara/h], were also decreased in Scap'^^^ mice. 

Suppression of the mRNAs of Wnt-pathway genes were 
confirmed by real time PCR analysis (Figure 7). Consistent with 
the RNA microarray findings, expression of several Wnt pathway 
genes were decreased in Scap'^'^ vs. control whole lung tissue at 
El 7.5 (Figure 7A). Interestingly, the response of Wnt signaling 
genes in the absence of epithelial Scap were dynamic with respect 
to age. Expression oiApc, Fzdl /2/ 4 and GskSb were slightiy higher 
in Scap'^'^ vs. control lung tissue at PNl (Figure 7B). We speculate 
that the suppression of Wnt pathway genes at El 7.5 was due to the 
direct effect of Scap deletion from alveolar epithelial cells, while the 
moderate increase in a subset at PN 1 may have resulted from the 
compensatory effect to the loss of SCAP in alveolar type II cells 
from adjacent cells such as lipofibroblasts, smooth muscle cells or 
endothelial cells [12]. 

To test if SCAP/SREBP signaling directly regulates Wnt 
pathway genes in pulmonary epithelial cells, we performed siRNA- 
mediated knockdown of Scap in an immortalized mouse lung 
epithelial cell line, MLE15 cells. This approach effectively reduced 
Scap mRNA expression to 10-20% of control levels (Figure 7C). 
Consistent with the microarray observation at El 7.5, siRNA 
knockdown of Scap reduced expression of all the Wnt pathway 
genes tested except Soxll and Lefl (Figure 7D). Soxll expression 
is undetectable in MLE15 cells, and results for Lefl expression 
were inconsistent using two independent Scap siRNA duplexes. 
Collectively, in vivo, in vitro and in sitico data support the concept 
that SCAP/SREBP signaling has stage-specific roles in the 
regulation of Wnt pathway genes during lung development and 
maturation. 

Suppression of Glucocorticoid Receptor Signaling by 
Epithelial Scap Deletion 

Glucocorticoid receptor (GR) signaling plays important roles in 
many physiological processes, including lung development and 
maturation [23]. Cellular responses to glucocorticoids are medi- 
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Figure 4. SREBP transcriptional network regulates lipogenic pathways. A) Clustering analysis of genes differentially expressed in Scap 
mice vs. control at El 7.5 and PN1 using a Self-organizing map. Genes from Clusters 5 and 7 were functionally enriched in lipid biosynthesis and 
metabolism. B) Upstream network analysis of Cluster 5 genes. Genes down-regulated in response to Scap deletion are shown as green nodes, genes 
that were unchanged but predicted to be involved in the process are shown in orange (activation) and blue (inhibition) respectively. C) Upstream 
network analysis of Cluster 7 genes. In addition to Srebfl/2, Cebpa, Foxa2 and Foxol were decreased in response to Scap deletion (shown as green 
nodes). Genes that were unaltered but predicted to be involved in the regulation of Cluster 7 genes, including Por and Pten (orange nodes indicate 
activation), Tgfbl, Nl<x2-1, Ifng, Pparg, Ctnnbl and Hifla (blue indicates inhibition), are shown. 
doi:10.1371/journal.pone.0091376.g004 



ated primarily througli binding to the intracellular GR (NR3C1) 
[24] and translocation of this complex to the nucleus where it 
interacts with multiple transcription factors including JUN, 
Nuclear Factor- KappaB, Octamer transcription factors OCTl 
and OCT2, CREB binding protein and STAT5 to positively and 
negatively modulate gene expression. Expression of multiple GR 
cofactors and target genes were altered in Scap^^^ lungs with 
respect to controls at E17.5, E18.5 and PNl (Table 5). The effect 
of SREBP on GR signaling was primarily seen at El 7.5, where 
most of the GR-associated genes were decreased in perinatal 
Scap'^^^ lungs. Similar to the alterations observed in Wnt signaling, 
a subset of genes displayed a moderate increase of expression at 
PNl, which may be due to compensation via adjacent cells. 
Glucocorticoid receptor {Nr3cl), three GR co-activators (Ncoal, 
Ncoa2, and NcoaS) and GR target genes [Tgfbl, Tgfh2 and SmadS) 
were significantly decreased in Scap'^'^ lungs at El 7., 5 and El 8.5. 



Nr3cl interacts with Ncoal and Ncoa2 to regulate transforming 
growth factor beta (Tgfhl) signaling via increased transcriptional 
activation of SMAD family member 3 [SmadS) [25-27]. Ptgs2 was 
one of four genes linked to GR signaling that was significantiy 
increased in perinatal Scap^^^ lungs. In the rat renal cortex, Ptgs2 
expression is suppressed by endogenous glucocorticoid levels 
through NrScl and inhibition of NrScl results in increased 
expression of Ptgs2 in the presence of dexamethasone [28-30]. 
These data indicate that epithelial Scap/Srebp signaling influences 
the expression of multiple GR signaling components in the lung in 
a time-dependent manner. 

miRNAs Responsive to Scap Deletion in Perinatal Lung 

MicroRNAs control many physiological and pathological 
processes through post-transcriptional regulation of gene expres- 
sion. To assess the role of the SREBP signaling cascade in the 
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Figure 5. Putative SREBP co-factors mediating perinatal lung lipid homeostasis. A SREBP related transcriptional network was constructed 
based on the promoter analysis of genes within Clusters 5 and 7. Each connection between a matrix family (white rectangle) and a target gene (blue 
oval) represents at least one predicted binding site for a given transcription factor on the promoter region of that target gene. 
doi:1 0.1 371/journal.pone.0091 376.g005 



regulation of miRNAs, genome-wide miRNA expression profiling 
was compared in Scap^'^ and control mice at El 8.5. We identified 
22 miRNAs that were differentially expressed in response to Scap 
deletion in epithelial cells (Table 6). Among these, 1 1 miRNAs 
were increased and 1 1 were decreased. Target genes were 
predicted for differentially expressed miRNAs using the miRNA 
Target Filter tool in Ingenuity Pathway Analysis. Target genes 
were further filtered if their differential expression patterns were 
negatively correlated to those of differentially expressed miRNAs 
in Scap^^"^ vs. control mice. Several miRs associated with the 
regulation of lipid metabolism were differentially expressed in 
Scap-^'" lungs at E18.5, including miR-122, miR-33, miR-29, 



miR-106 and miR-335 (Table 6). These data suggest that, in 
addition to the transcriptional regulation of lipid homeostasis by 
SREBPs, altered miRNAs in response to Scap/Srebp deletion is 
likely to further modulate lipogenic gene expression at the 
posttranscriptional level to regulate lipid homeostasis. 

Several miRNAs known to influence lung development, 
maturation and disease were also differentially expressed in 
Scap'^''^ vs. control lungs including miR-l, miR-124, miR-29 
and miR-l 33. Previous data demonstrated a link between miR-l 
and glucocorticoid signaling. Treatment with dexamethasone and 
myostatin induced miR-l expression via nuclear translocation of 
NrScl [31]. miR-l is known to be regulated by Srebfl and Cehpa, 
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Table 3. Potential upstream regulators of SCAP/SREBP pathway (ranked by P-value). 



E17.5 Only 




El 7.5 Insig"" 




El 7.5 Scap*" 




El 8.5 Scap'^"^ 




PN1 Scap'"* 




Gene 


P value 


Gene 


P value 


Gene 


P value 


Gene 


P value 


Gene 


P value 


F0X03 


u.uuuDyy 


SREBF2 


2 44E-09 


WTl 


5 33E-06 


SREBF2 


3 26E-06 


SREBF2 


1 02E-1 2 


IRF7 


n nnnQQ7 
u.uuuyy/ 




2 07E-08 


SPDEF 


7 2E-06 


F0X04 


4 1 7E-06 


SREBFl 


3 51 E-1 0 


E2F4 


0.00373 


SIRT2 


5.33E-06 


PROXl 


2.32E-05 


SREBFl 


1.85E-05 


SIRT2 


3.59E-10 


CTNNBl 


0.00663 


TP53 


1 .76E-05 


KLF2 


7.18E-05 


NR1H3 


0.000022 


F0X04 


4.6E-07 


P0U4F1 


0.00881 


PPARA 


0.000106 


HHEX 


9.65E-05 


SIRT2 


3.64E-05 


PPARA 


2.49E-05 


GUI 


0.00989 


WTl 


0.000106 


NFE2L2 


0.000435 


NR1H2 


0.000452 


NR1H2 


7.24E-05 


GATA3 


0.0231 


STAT3 


0.00018 


NOTCH 1 


0.000699 


STAT6 


0.000464 


PPARGCl B 


0.000173 


SPDEF 


7.2E-06 


Smad2/3/4 


0.000599 


HTT 


0.000829 


PPARGCIB 


0.00071 


NR1H3 


0.000189 






TOBl 


0.000724 


NFkB 


0.000905 


PPARA 


0.00305 


PPARG 


0.00088 






SREBFl 


0.000742 


IRF7 


0.000997 


ETSl 


0.00485 


ESRRA 


0.0017 
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both of which are important regulators of lipogenesis in alveolar 
type II cells [32]. miR29 is associated with the pathogenesis of 
pulmonary fibrosis [33] and regulates pro-inflammatory cytokine 
secretion by targeting Insigl, Srehfl and Lpl [34,35]. miR-133, 
regulated by SREBP 1 [32], is associated with inflammatory lung 
disease [36]. Our data support the concept that a group of 
microRNA involved in lipid metabolism and lung development/ 
maturation were dysregulated in response to Scap deletion in 
epithelial cells, potentially contributing to the crosstalk between 
SREBP signaling and pathways critical for lung development and 



maturation (e.g., glucocorticoid signaling pathway), connecting the 
"lipid" and "lung" related transcripts within the SREBP network. 

Discussion 

We sought to identify the role of epithelial SREBP, a known 
regulator of cellular lipid homeostasis in multiple tissues, during a 
critical time period of perinatal lung development and maturation 
from El 7.5 to PNl. Our data show that SCAP/INSIG/SREBP 
signaling primarily influences epithelial development, cell death 



Table 4. Potential upstream regulators of SCAP/SREBP pathway (ranked by Z-score) 



Inslg'^" 
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Gene 


State 


Z-score 
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State 


Z-score 
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4.1 
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4.464 


SREBFl 


Activated 


3.999 


PPARGCl A 
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3.683 


CDKN2A 
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3.805 


SMAD7 


Activated 


3.407 


STATl 
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3.457 


NRFl 


Activated 


3.283 


SREBF2 
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3.375 


SPDEF 


Activated 


3.099 


IRFl 


Activated 


3.372 


HICl 


Activated 


2.985 


SPDEF 
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3.046 


HNF4A 
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2.937 


AHR 


Activated 


3.022 


ESRRA 


Activated 


2.852 


CEBPA 
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2.972 


CEBPA 


Activated 


2.244 


TRIM24 
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-3.065 


SREBFl 


Inhibited 


-5.146 


NC0A3 


Inhibited 


-2.932 


SREBF2 


Inhibited 


-4.227 


HOXAIO 


Inhibited 


-2.797 


PPARGCIB 


Inhibited 


-3.614 


FOXMl 


Inhibited 


-2.705 


NRIPl 


Inhibited 


-3.187 


E2f 


Inhibited 


-2.508 


NFYA 


Inhibited 


-3.172 


BCL6 


Inhibited 


-2.497 


E2f 


Inhibited 


-2.795 


TCF4 


Inhibited 


-2.312 


Creb 


Inhibited 


-2.643 


NRIPl 


Inhibited 


-2.262 


ESRl 


Inhibited 


-2.621 


MYC 


Inhibited 


-2.258 


Nrlh 


Inhibited 


-2.59 


MEF2D 


Inhibited 


-2.246 


GATA3 


Inhibited 


-2.388 
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Figure 6. Epithelial Scap deletion in the perinatal lung suppresses Wnt/p-catenin signaling. (A) A schematic representation of mRNAs 
changed in the canonical Wnt/)3-catenin pathway in Scap^^'^ perinatal lung (El 7.5) is shown. Green nodes indicate mRNAs that were decreased in 
Scap'^^^ lungs with respect to control lungs. Wnt pathway genes that were not changed in response to Scap deletion are depicted as empty nodes. 
(B) Heatmap demonstrates decreased expression of genes encoding a number of key components in the Wnt/p-catenin pathway in Scap^^'' vs. 
control mice, n = 3/group. 
doi:10.1371/journal.pone.0091376.g006 



and cell proliferation at El 7.5, whUe primarily influencing 
surfactant physiology, lipid/sterol synthesis, and phospholipid 
transport at PNl. In addition to regulating the lipogenesis 
program that was dominant at the later stage of lung development, 
this study demonstrates cross-talk between SREBP signaling and 
two other signaling pathways important for lung development and 
maturation, Wnt/ P-catenin and glucocorticoid receptor signaling. 
Epithelial SCAP/INSIG/SREBP signaling regulates expression of 
miRNAs that control lipid homeostasis and glucocorticoid 
signaling at the post-transcriptional level. Taken together, our 
analysis demonstrates that SCAP/INSIG/SREBP regulates peri- 
natal lung lipogenesis and maturation through its interactions with 
diverse regulatory partners. 

SatPC pool sizes were significantiy decreased in Scap'^ lungs at 
E18.5 and PNl (Figure 2A). SatPC production and secretion by 
alveolar type II cells increases prior to birth to facilitate the 
transition to air breathing. Synthesis of SatPC occurs by direct de 
novo synthesis (Kennedy pathway) and the remodeling of existing 
monosaturated PC species in alveolar type II cells (Lands cycle, 
also known as the remodeling pathway); the relative contribution 
of these pathways to the total SatPC pool is approximately equal in 
vitro and in vivo Pcytla is the rate-limiting enzyme for the de novo 
synthesis pathway and is critically required for SatPC production 



and lung function at birth [37] . Two genes in the de novo synthesis 
pathway [Pcytla and Pld4) and two in the remodeling pathway 
[Lpcatl and 3) were significantly decreased in perinatal Scap^^^ 
lungs, all four of which belong to the lip id-enriched cluster. Cluster 
.5 (Figure 4B, Table SI). Lpcatl and LpcatS encode lysopho- 
spholipid acyltransferases that, in conjunction with phospholipas- 
es, generate SatPC via the remodeling pathway. Lpcatl primarily 
reacylates lysoPC with a saturated fatty acid donor and is required 
for pulmonary SatPC production in fetal mice [38]. LpcatS also 
demonstrates lysoPC acyltransferase activity, is expressed in the 
lung and at high levels in the liver where it drives phospholipid 
remodeling [39,40]. Pcytla, Lpcatl and Pld4 have putative Srebp 
binding sites in their proximal promoters that are conserved 
between mouse, rat and human, suggesting that they are likely 
direct targets of the SCAP/SREBP/INSIG pathway (data not 
shown). Collectively, in addition to driving the initial phases of 
lipogenesis in the lung (discussed below), these data are consistent 
with the concept that Scap also drives the final steps of SatPC 
synthesis during late pulmonary maturation, in part, through 
induction oi Pcytla, Lpcatl, LpcatS and Pld4 expression. 

We used a combination of profile clustering, promoter mining 
and literature mining approaches to identify potential upstream 
regulators and downstream targets of SCAP. Through this 
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Figure 7. 5ca/3 modulates expression of Wnt pathway components. A and B) qPCR analysis demonstrates altered expression of Wnt signaling 
genes in Scap'''^ El 7.5 (A) and PN1 (B) whole lung tissue compared to littermate controls. Results were normalized to 18 s RNA, n = 3-4/group, *p< 
0.05 vs. control at the respective time point. C and D) qPCR analysis demonstrates decreased expression of Scap (C) and of Wnt pathway (D) genes in 
IV1LE15 cells following transfection with negative control (neg con) and two distinct Scap siRNAs (sil and si2). Results were normalized to 18 s, n = 6 
samples/group representing 3 independent experiments, *p<0.05 vs negative control. ND = non-detectable. 
doi:10.1371/journal.pone.0091376.g007 



analysis, difTerentially expressed genes in Scap'^'^ vs. control mice 
were placed into 8 distinct clusters based on expression similarity. 
Genes in clusters 5 and 7 were both functionally involved in the 
lipid synthesis, metabolism and transport. Additionally, cluster 7 
genes also play important roles in lung development and function 
including Sftpb, AbcaS, FgflS and Pdgfa. Mutations in the genes 
encoding Sftph and AbcaS are associated with respiratory distress 
and interstitial lung disease in the pediatric population (SFTPB; 
Online Mendelian Inheritance in Man [OMIM] number 178640, 
ABCAS; OMIM number 601615) and mutations in the human 
SFTPAl gene (c.655C>T (rs4253527)) increase the risk of 
idiopathic pulmonary fibrosis (http:/ /www.ncbi.nlm.nih.gov/ 
clinvar/RCVOOOO 14093.1/). FgflS is required for late stage 
embryonic lung alveolar development as FgflS deficiency in mice 
resulted in reduced alveolar spaces associated with thickened 
interstitial mesenchymal compartments [41]. Pdgfa is critical for 
alveolar myofibroblast ontogeny. Postnatally surviving Pdgfa(— /—) 
mice develop emphysema secondary to the failure of alveolar 
septation [42]. Other notable genes in cluster 7 include Hesl, 
Hyall, Clic5, Cdol, Hck, and Twist2, the deletion/mutation of 
which all lead to abnormal pulmonary alveolar morphology 
(information from http://www.informatics.jax.org/phenotypes. 
shtml MP:0002132). 

Potential upstream regulator analysis of cluster 5 (C5) and 
cluster 7 (C7) genes identified SREBP as a common key regulator 
for both C5 and C7 genes; however, other potential key regulatory 
partners in C5 and C7 were distinct. As shown in Figure 4B, 
Hnf4a, Pax3, and J{fe2l2 were predicted as major regulators 
working in concert with SCAP/INSIG/SREBP signaling to 
control expression of C5 genes. Hnf4 (Hepatocyte nuclear factor 
4-alpha) was ranked as the 6th important key regulatory hub of C5 
based on over-representative p-values (SREBP2, SCAP, SREBPl, 
INSIGl, INSIG2 were listed 1-5) and as the top TF in the 
network based on the target gene numbers in C5 (regulating 15% 



of C5 genes, data not shown). HNF4 regulates a number of genes 
involved in fatty acid synthesis {Acly, Asahl, Midlipl, Scdl and 
Tp53); phospholipid synthesis (Acp6, Pcytla, Lypla2, Pitpnh and 
Scdl); and lipid transport {Slc25al, Slcl7a5, Rab7a, and Asahl). 
Multiple lines of evidence suggest that SREBP 1/2 directly binds 
and regulates Hnf4 expression in vivo and in vitro [43-46]. Nfe2l2 
(Nuclear factor erythroid 2-related factor 2) and Pax3 (paired box 
gene 3) influence Srebfl expression in mouse liver and human 
meduUoblastoma cell line [47,48]. On the other hand, Por, Pparg, 
Cehpa, Nkx2-1, Tgfbl, Ctnnbl, Foxol, Foxa2, and Ifiig were predicted 
as central hubs working together with SREBP to regulate gene 
expression in cluster 7 (Figure 4C). Srebfl, Cehpa and Pparg are 
known transcription factors that regulate lipogenesis in various 
tissues. CEBPA directly binds to the promoter of Srebp-lc and 
positively regulates SREBP Ic expression during adipogenesis [49]. 
Activation or deletion of SREBP Ic lead to a corresponding 
increase or decrease in CEBPA expression in hepatocytes and 
adipocytes [49-51]. Similarly, SREBPlc and PPARG transacti- 
vate each other's expression in hepatic cells and adipocytes [52- 
55], therefore SREBP, CEBPA and PPARG share the same 
transcriptional network controlling lipogenesis and metabolism. 
Foxol is a master regulator of energy metabolism in multiple tissues 
[56]. In the liver, Cehpa forms an interaction complex with Foxol to 
suppress Foxol expression [57,58] and Foxol inhibits SREBP-lc 
expression and other lipogenic genes [59]. Previous work has 
shown that NKX2-1 and SREBP co-regulated Abca3 expression 
through two distinct cis-acting regulatory cassettes, a 'lung 
regulatory unit' (TTF-1, GATA, FOXA2, CEBPA, etc.) and a 
'lipid regulatory unit' (SREBP, SPl) present in AbcaS promoter 
[60] . It is likely that genes sharing both 'lung' and 'lipid' properties 
in cluster 7 are regulated through similar mechanisms. 

One unexpected finding from this study was suppression of 
multiple W nt pathway components in the Scap'^ lungs, including 
ligands [Wnt 5/5a/7d), receptors (Fzd 11214) and effectors [Gsk3h, 
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Table 5. Glucocorticoid receptor signaling 


was suppressed 


In Scap^'^ mice 


in perinatal lung at El 7.5. 




Symbol 


E17.5 


E18.5 


PN1 


RefseqID 


Foxo3 


052 


0.79 


0.85 


NM_019740 


Polr2a 


0.54 


0.81 


1.27 


NM_009089 


Map3kl 


058 


0.77 


1.03 


NM_011945 


Nos2 


0.60 


0.79 


0.68 


NM_0 10927 


Pik3r4 


0.60 


0.87 


0.95 


NM_0C1081309 


Smad3 


0.61 


0.80 


0.93 


NM_016769 


Tgfb2 


0.62 


0.77 


0.85 


NM_009367 


Pou2fl 


0.63 


0.77 


1.52 


NM_011137 


Polr2b 


0.66 


0.90 


0.73 


NM_1 53798 


Sntiarca2 


0.66 


0.90 


1.01 


NM_011416 


Smarca4 


0.66 


0.77 


1.30 


NM_001 174078 


Crebbp 


0.66 


0.79 


1.82 


NM_001025432 


Nripl 


0.67 


0.87 


1.19 


NM_1 73440 


Jun 


0.67 


0.82 


0.89 


NM_010591 


Atm 


0.67 


0.90 


1.19 


NM_007499 


Nfatc3 


0.67 


0.83 


0.75 


NM_010901 


Ncoal 


0.68 


0.85 


1.24 


NM_0 10881 


Sosl 


0.68 


0.90 


0.91 


NM_009231 


Ncor2 


0.68 


0.78 


1.41 


NM_0 11424 


Nfkbl 


0.69 


0.92 


0.84 


NM_008689 


Pik3rl 


0.69 


0.84 


1.38 


NM_0C1 077495 


Ncoa3 


0.70 


0.88 


1.10 


NM_008679 


Cdk7 


0.71 


0.90 


058 


NM_009874 


Nr3cl 


0.72 


0.87 


0.97 


NM_008173 


Ncoa2 


0.76 


0.86 


1.43 


NM_008678 


Mapkl4 


0.80 


0.89 


0.87 


NM_011951 


Mapk3 


0.81 


0.89 


0.68 


NM_011952 


Cebpa 


0.82 


0.91 


0.48 


NM_0C7678 


StatSb 


0.84 


0.88 


0.71 


NM_0C1 113563 


MapklO 


0.91 


1.13 


0.68 


NM_009158 


Pik3c2g 


1.15 


1.32 


1.83 


NM_207683 


Polr2k 


1.40 


1.84 


1.30 


NM_001 039368 


Ptgs2 


1.98 


1.72 


0.50 


NM_011198 


doi:l 0.1 371 /journal.pone.0091 376.t005 



Ctnnbl) (Figure 6, Figure 7A and 7B). Results initially identified by 
microarray analysis were confirmed at the niRNA level in MLE 1 5 
cells in which endogenous Scap was inhibited, demonstrating a 
direct efiect of Scap in the regulation of specific Wnt signaling 
pathway components (Figure 7D). Wnt signaling is critically 
required for pulmonary epithelial and vascular development [61- 
64] and has been implicated in chronic lung diseases including 
pulmonary fibrosis [65-67] and chronic obstructive pulmonary 
disease [68,69] . Defects in lung morphology were not observed in 
Scap'^'^ lungs in the present study (data not shown). One possible 
explanation is that the degree of Wnt suppression in Scap'^''^ lungs 
was insufficient to phenocopy mouse models of epithehal Wnt loss 
[70,71]. It is also possible that residual Scap expression, due to 
incomplete deletion of Scap from epithelial cells in our model, was 
sufficient to maintain Wnt signaling adequate to support normal 
pattering of the lung. 



In addition to its role in organogenesis, the Wnt signaling 
pathway has emerged as a key regulator of lipogenesis in mice and 
humans where it acts as a suppressor of adipocyte maturation. 
Three Wnt ligands have been implicated in suppression of 
adipogenesis in mice including WNT 10b, WNT5a and WNT5b 
[72]. Enhanced Wnt signaling blocks adipogenesis via inhibition of 
the central transcriptional regulators CEBPA and PPARg [73]. 
The biological significance of Wnt suppression in Scap'^ lungs is 
currently unknown. However, we postulate that suppression of 
Wnt signaling in Scap'^'^ alveolar type II cells serves as a 
compensatory mechanism to increase surfactant lipogenesis in the 
absence of SCAP/INSIG/SREBP signaling. The compensatory 
regulation of lipogenesis in the absence of SCAP in epithelial cells 
is consistent with previous findings in the adult lung wherein 
phospholipid synthesis was decreased in type II cells of adult 
Scap'^'^ mice; storage, synthesis, and transfer of lipids by lung 
lipofibroblasts were increased [12]. Similarly, SCAP deletion in 
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Table 6. Differentially expressed miRNA in El 8.5 Scap vs. control mouse lungs. 



Symbol 


Fold Change 


Target Count 


Function and Reference 


miR-122 


2.5 


17 


l<ey regulator of cholesterol and fatty-acid metabolism (Esau, 2006); regulate lipogenic gene expression 
including SREBF-lc, Fasn, Hmgcsl, Sqle and Acaca (Krutzfeldt, 2005; lliopoulos, 2010} 


miR-196a 


1.79 


23 


inducible brown adipogenesis by suppresse the expression of the white-fat gene Hoxc8 (Mori, 2012) 


miR-602 


1.75 


10 




miR-449c 


1.72 


51 


highly expressed in lung shortly before and after birth, induce epithelia differentiation by repressing the 
Notch pathway (Bao, 2012, Lize, 2010) 


miR-1202 


1.64 


27 




mir-4784 


1.64 


4 




miR-647 


1.56 


32 




miR-29a 


1.54 


21 


pathogenesis of pulmonary fibrosis (Gushing, 2010), regulates pro-inflammatory cytokine secretion by 
targeting Inslgl, Srebfl and Lpl (Chen, 201 1,Hulsmans 2011) 


miR-192 


1.54 


1 




miR-15a 


1.52 


8 


pre-adlpocyte differentiation, targeting DLKl (McGregor, 2011) 


miR-375 


1.28 


7 


alveolar epithelial cell differentiation by inhibiting Wnt/p-catenin pathway (Wang, 2012) 


miR-103a 


-1.3 


15 


regulation of lipid metabolism (Wilfred, 2007) 


miR-33 


-1.35 


8 


key post-transcriptlonal regulator of cellular cholesterol and lipid homeostasis, hosted within the introns 
of Srebf2 and Srebfl (Rayner, 2010; Marquart, 2010; Najafl-Shoushtari, 2010} 


miR-335 


-1.37 


2 


up-regulation of microRNA-335 is associated with lipid metabolism and obesity (NakanlshI, 2009) 


miR-363 


-1.54 


22 




miR-185-5p 


-1.57 


17 


regulates adipocyte differentiation, up-regulated in mature adipocytes while down-regulated In obese 

IMcll ^Vw'ILcga, 1 


miR-376b 


-1.58 


26 




miR-204 


-1.62 


29 


regulated by STAT3, down-regulated in pulmonary arterial hypertenslon(Gourboulin, 2011) 


miR-1 


-1.66 


29 


regulate by SREBF1 and Gebpa (Granjon, 2009) 


miR-290-5p 


-1.8 


12 




miR-1 24-3p 


-1.86 


43 


down-regulated during fetal lung development. Inhibits fetal alveolar epithelial cell maturation 


miR-1 33b 


-1.91 


25 


regulate by SREBFl (Granjon, 2009}, inflammatory lung disease(Oglesby, 2010); Insulin-dependent 
diabetes mellltus (Li, 2009) 



doi:10.1371/journal.pone.0091376.t006 



the liver resulted in a compensatory increase in fatty acid synthesis 
in adipocytes [74]. 

MicroRNAs have emerged as key regulators of various cellular 
metabolic pathways, including lipogenesis. In this study we 
identified multiple miRs that were differentially expressed in 
Scap^^^ mice. miR122 and miR33 are known to directly regulate 
lipid metabolism and lipid homeostasis in the liver [75,76]. miR- 
33 was identified as a key post-transcriptional regulator of cellular 
cholesterol homeostasis and is present in several animal species, 
including humans. Humans possess two family members called 
mir-33a and mir-33b, which are located within intron 16 oi Srebf2 
and intron 17 ot Srebfl, respectively. Both mir-33a and mir-33b are 
co-transcribed with their host genes {Srebfl/ 2) and reciprocally 
regulate genes involved in cholesterol and fatty acid metabolism 
[77-79]. Probes targeting multiple species of miR-33 were 
decreased in Scap^'^ vs. control mice. miR-122 is highly conserved 
from human to frogs and most abundantly expressed in liver. 
Blockade of miR-122 expression decreased Srebf-lc, Fasn, Hmgcsl, 
Sqk and Acaca expression, resulting in the reduction of plasma 
cholesterol and triglyceride levels in both rodents and primates 
[80,81]. The fact that miR-122 expression correlated with Srebfl 
and other lipid metabolism related genes indicated that these genes 
are not direct targets of mLR-122. The mechanism by which 
miR122 activates these genes and lipid metabolism is currently 
unknown. miR122 was induced 2.5 fold in Scap^''^ vs. control 
mice, perhaps serving as a compensatory response to modulate 



SREBP and other lipogenic genes. Other miRNAs involved in 
regulating lipid homeostasis including miR-29 [82], miR-1 06 [83], 
and miR-335 [84] were differentially expressed in Scap^'^ vs. 
control mice. Our data suggest that in addition to the transcrip- 
tional regulation of lipid homeostasis by SREBPs, miRNAs may 
play important roles at the posttranscriptional level to influence 
lipogenic gene expression and maintain lipid homeostasis. 

In summary, our previous network model identified the 
important role of SREBP as a regulatory hub controlling a lung 
lipid transcriptional network [14]. The present study was designed 
to further define the dynamic roles of SREBP during perinatal 
lung development and maturation. SREBP signaling influences 
genes controlling epithelial development, cell death and cell 
proliferation at El 7.5 and those involved in surfactant physiology, 
lipid/sterol synthesis, and phospholipid transport at a later 
gestational age (PNl). Our study supports the concept that 
SREBP signaling interacts with other signaling pathways impor- 
tant for lung development and maturation, including the Wnt/fS- 
catenin pathway. SREBP regulates perinatal lung lipid homeosta- 
sis and lung maturation through multiple mechanisms including its 
interactions with diverse regulatory partners to regulate target 
genes and microRNAs that are known to regulate lipid metabolism 
at the posttranscriptional level. 
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